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ABSTRACT 



The power spectral density of a pulse ~po sition modulated (PPM) 
waveform is determined in part by certain statistics of the modulating 
process. In this paper a mathematical model is developed to determine 
the effects of three different modulating processes upon the power spec- 
tral density. The processes considered are the Gaussian, the exponen- 
tial, and the uniform processes, with statistical independence assumed 
between the sample values that modulate the PPM waveform. The 
Gaussian process is also developed for the case in which statistical in- 
dependence may not be assumed; that is, a correlation exists between 
the sample values of the modulating process. The model will also 
determine the power spectral density of pulse -amplitude modulated 
(PAM) waveforms, or for a combination of PAM and PPM, 

The mathematical model is reduced to a computer model and 
computed results are displayed. The design of a pulse -po sition modu- 
lator is shown, and experimental measurements of power spectra are 
compared with the theoretically determined spectra. The theoretical 
determination of power spectral density and the actual measurement of 
spectra of finite-time samples are discussed. 
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the writer during his Industrial Tour at Sylvania’s Electronic Defense 
Laboratories (EDL) in Mountain View, California. 

The writer gratefully acknowledges the helpful assistance and 
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Introduction 



The objective of this report is to relate a theoretical approach 
for the determination of power spectral density, to show computations of 
power spectral densities performed for a particular class of signals, and 
to compare experimental measurements of the power spectral densities 
to that predicted by theory. 

The power spectral density of a pulse -po sition modulated (PPM) 
waveform typically consists of a continuous component and discrete fre- 
quency components. The relative amount of the total power contained in 
each of the component types is determined by the modulation process. 

The distribution of power in the spectrum is determined by the modula- 
tion process, the pulse shape, and the average pulse repetition rate. 

In this study, a mathematical model is developed which determines 
the power spectral density of PPM waveforms. The mathematical model 
is then reduced to a digital computer program. The computer program 
permits the computation of the power spectral density of a PPM wave- 
form that is modulated by a uniform process, a Gaussian process, or an 
exponential process. The modulation process that results in a statistical 
correlation between pulses within the pulse train with respect to their 
displacements (in time) is developed for the Gaussian process only. Sta- 
tistical independence between pulses is assumed for the other processes 
considered. 

The pulse shape is assumed to be rectangular for the computer 
model. The pulse width and the average pulse repetition rate are variable 
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parameters of the program. The pulse amplitude is also a variable param- 
eter, thereby allowing computations to be performed for combined pulse- 
position and pul se -amplitude modulated waveforms. By setting the pulse 
position modulation equal to zero, the computation of the power spectral 
density of a pulse waveform that is only pulse-amplitude modulated may 
also be accomplished with the computer model. 

A discussion of the use of finite-time interval measurements, 
using a periodogram technique to give estimates of the power spectral 
density of a process, is included in this report. The periodogram tech- 
nique is used to determine the power spectral density of a PPM waveform. 
Since the characteristics of the experimental PPM waveform are also 
used as an input to the computer model, a comparison of the theoretical 
and experimental results is included. 
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2 . 



A Discussion of Spectral Analysis (General) 



The purpose of a communication system is to convey information 
in some form from one place to another. The system would normally 
include a transmitter, a transmission medium, and a receiver which 
recreates a reasonable representation of the original message. The 
message creating the transmitted signal will most likely be one of a 
class of stochastic phenomena; that is, the message is a random func- 
tion of time. Also, the noise accompanying the signal (introduced by 
the communication system) will normally be a stationary Gaussian ran- 
dom process. The analysis of such signals is difficult unless one re- 
sorts to the use of the statistical properties of the signals. The power 
spectral density furnishes us such an analysis technique, as it is rela- 
ted to the statistical parameters of the signal. In particular, the power 
spectral density is related to the covariance (or autocorrelation) function. 

In the development of spectral analysis techniques, one normally 
considers the random signal to be a stationary process. A stationary 
process is one in which the statistical parameters of the signal are the 
same everywhere, independent of which long interval over which they 
are measured. A similar concept is to consider a process as an en- 
semble of signals (an ensemble is a set of similar things). If the statis- 
tical averages of the ensemble measured at one instant of time are the 
same for any other instant of time, then the process is also said to be 
stationary in time. 

If a single process of the ensemble of signals has statistics over 
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a long period of time that are the same as the statistics of the ensemble 
at a given time, then the process is said to be "ergodic". An ergodic 
process will always be stationary, but a stationary process does not 
necessarily imply ergodicity [6] . 

In this paper all processes will be assumed stationary unless 
otherwise stated. 

To arrive at the statistical parameter used to obtain the power 
spectral density, consider the random ergodic process y(t). A typical 
member of the ensemble is y n (t). The ensemble averaged correlation 
function (covariance function K(T))is defined as the second joint moment, 
or 

K(T) = E {V(t + T)} ’ 

where E denotes the expected value taken over the ensemble. ^ The auto ~ 
correlation function R n (T) of a typical member of the ensemble is 

_n. . lim 1 pT n, . n, t . , 

R (T) = T~co If J_ T y (t)y (t+ T)dt 

In general, 

K(t) = E ^R n ( t ) j- , 

where the E-j^R n (T)j- here is the autocorrelation function averaged over the 

The statistical definition of covariance would require that the 
square of the mean (tiy^) be subtracted from the above expression. It was 
decided, however, to include the mean in this investigation, thus permit” 
ting a simplified development that will handle signals that possess a non- 
zero mean. 
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entire (infinite) ensemble. This average should normally be performed 
before the time average of R(T). 



However, for the ergodic case 
R n (T) = K(T) (probability 1) . 



The mathematics relating the power spectral density to the auto co rrela- 
tion (or covariance) function were first developed by Wiener and 
Khintchine. The Wiener -Khintchine Theorem states that the power 
spectral density S(f) is the Fourier transform of the covariance func- 
tion, or 



S(f) 



poo 

= z 1 
“OO 



K(T) e “^ U,T dT > f ^ o 



Note that in the above expression S(f) is defined for positive frequencies 
only, hence the factor of two in the equation. This is obviously the more 
practical case, and is permis sable since S(f) is an even function (as is 
R(T)). One can obtain the autocorrelation function by taking the inverse 
Fourier transform of the power spectral density (again assuming ergo- 
dicity). Thus 



i r°° 

R(T) = •— S(f)cos uoTdf 
o 



(probability 1) 



and 



pOO 

S(f) = 4 R(T)cos uuTdT 
o 



(probability 1) 



It should be mentioned that while a signal has a unique power spectral 
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density, the converse is not true. All phase information in the signal 



is lost in S(f), thus the original signal could not in general be recov- 
ered from its power spectral density. Furthermore, a particular S(f) 
theoretically represents the power spectral density for a multitude of 
signals. 
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3. 



Finite-Time Interval Spectral Analysis 



The developments of the previous section followed the typical 
theoretical approach to signal analysis in that infinite time intervals 
were assumed. In practical analysis situations, such intervals are of 
course always finite. For some processes, a finite investigation in- 
terval can be considered adequate if the interval is large with respect 
to the fluctuation time of the process under investigation. Since finite 
time measurements of spectra of PPM waveforms will be compared 
with the results from the mathematical model, a discussion of possible 
sources of error due to analysis of the finite-time interval observation 
of the signal is included. 

The problem of finite time interval analysis of a stationary ran- 
dom process can be stated quite simply as follows. A sample function 
of the random process y(t) is observed during an interval t - -T/2 to 
t = T/2. Will this observation contain enough of the characteristics of 
the process to give a reasonable estimate of the power spectral density 
of the process? Note that the random process under consideration is 
stationary. The stationarity of the signal becomes more meaningful 
(and necessary) when considering finite-time analysis. 



the ensemble y(t), except now let y(t) vanish outside the interval (-T/2, 
T/2). An amplitude spectrum for a member of the ensemble is 



To develop the theory for finite -time analysis, again consider 
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Now the finite-time power spectral density of the ensemble is 

s(f) T = E {|- Y n (f) T Y n (f) T * } = ~ E { | Y n (f) T l 2 } , 

❖ 

where denotes the complex conjugate. The true power spectral den~ 
sity S(f) is defined in terms of the ensemble average over an infinite 
time interval (and is equivalent to the definition on page 5). Specifically , 

S(f) = T-co S(f) T = ^ E{|Y n (f) T | 2 } 

where the ensemble average is carried out before the limit is taken. 

Again, consider the sample function of a single random process. 
The finite time transform is equivalent to Y n (f)^, of the preceding para- 
graph, and will now be used as an approximation to the Fourier trans- 
form of a single stationary process. A finite power spectral density of 
the single random signal may now be defined. 

S n (f,T) = £ | Y n (f) T | 2 

2, 3 

S(f, T) is called a n periodogram !I , * and has received considerable 
attention due to the fact that it is one type of practical measurement 
that can be made on real processes (signals). Observe that the power 

Middleton, D. An Introduction to Statistical Communication 
Theory. McGraw-Hill, i960, 143 

Abramson, N. and Farison, J. On Statistical Communication 
Theory. Technical Report No. 2005-1, SEL-62-078, Stanford Univer- 
sity, August 1962* 30-33 
3 

Grenander, U. and Rosenblatt, M. Statistical Analysis of 
Stationary Time Series. John Wiley & Sons, 1957, 91-94 
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spectral density for the time limited ensemble is the expected value of 
the periodogram over the entire ensemble, or 

S(f) T = E{s n (f s T)} 

and 

S<« = E{s n (f,T>} 

Thus, one measure of the reliability of the periodogram is the fact that 

as T-*oo, the mean value of the periodogram does equal the power 

spectral density. However, as shown by Abramson and Farison, for 

some very common random processes, the variance of S(f, T) does not 

2 

go to zero as T -* oo. The fact that the periodogram (being a random 
variable) has the desired mean value but may vary extensively about the 
mean for certain processes has caused some skepticism about its use 
as a spectral density estimator [ 5] . 

It has been shown that the periodogram is not a consistent esti- 
mate of the power spectral density [10]* However, the periodogram 
may be used to give consistent estimates of the power spectral density 
at specific frequencies f . This is accomplished by the use of weight- 
ing functions, or n spectral windows’ 1 about the frequency f.. A possible 
method then of obtaining a consistent estimator of the power spectral 
density would be to have many adjacent windows, and at the end of the 
analysis time T observe the response of all windows simultaneously. 

A smoothing between the point estimates observed would naturally be 
desirable. 
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Periodogram calculators fall into three general classes; that is, 
the filtering, computing and smoothing may be performed in three dis- 
tinctly different ways. These are 

1. Parallel filter banks, 

2. Sweeping analyzers, 

3. The coherent memory delay integrator. 

The parallel filter bank technique is, as the name implies, a bank 
of narrow-band filters, each with a common input (the signal to be ana- 
lyzed). At the output of each filter is an energy detector, which is 
periodically sensed by a sweep circuit that is connected to the vertical 
deflection plates of an oscilloscope. The parallel filter bank technique 
is depicted in block diagram form in Figure 3-1. 




Figure 3-1. Parallel Filter Bank' Periodogram Calculator. 



The number of filters (N) may vary with the frequency range (analyzer 
bandwidth) and the resolution required. The bandwidth of a filter is 
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approximately the inverse of the rise time. Obviously, the bandwidth 
of the filters determines the resolution of the device. The ’’analysis 
time”, or Period T, using this method is approximately the rise time 
of the filter. From this, one may observe that the resolution and the 
analysis time are essentially the inverse of each other. Consequently, 
if one designs an analyzer for fine resolution, then the signal must be 
observed for a longer period, and the converse is also true. This im- 
portant relationship holds for periodogram calculators in general. 

The sweeping analyzer is the more common general laboratory 
type of device. With this device, the input signal is mixed with a sweep- 
ing local oscillator, the output of which is filtered through a single 
’’spectral window”. The signal is then detected and displayed on an os- 
cilloscope. A block diagram for the sweeping analyzer is shown in 
Figure 3-2. 




Figure- 3-2. Sweeping Analyzer Periodogram Calculator. 

A disadvantage of this type of analyzer is the adverse effect of the tran 
sient condition at the input to the filter while the local oscillator is 
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sweeping. If the sweep rate is too rapid, the filter may continue to ring 
due to this transient effect, even though the signal is no longer present 
at the input to the filter, and the resolution is no longer the bandwidth 
of the filter. 

The third type of periodogram calculator to be discussed is the 
coherent memory delay integrator. This device is somewhat sensitive 
in its operation, but the concept is unique and merits consideration. 

The block diagram of the device is shown in Figure 3-3. 



Readout after 




Output 



Figure 3-3. Coherent Memory Delay Integrator 



The local oscillator used for the phase advance signal must be precisely 

matched to the delay of the delay line. Furthermore, the loop gain must 

be unity. With these optimum conditions assumed in being, the response 

of the analyzer to a single frequency is a main lobe, with side lobe pat- 

2 

terns very much like a (sinx/x) curve. The side lobe patterns may 
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be removed in some instances by a process of selective "weighting 11 , 
in which the input signal is time weighted during the analysis interval. 
As am example, a cosine weighting function is accomplished by ampli- 
tude modulating the input signal. 

The output is produced by a linear addition of a large number of 
sinusoids differing from each other in frequency by small increments. 
Phase coincidence of all signals present in the storage loop occurs once 
(for a given frequency) during the observation interval. The time of 
this phase coincidence is linearly related to the input frequency, thus a 
readout of the response after the proper number of loop recirculations 
yields the desired periodogram. The one operational analyzer of this 
type with which this writer is familiar is known under the trade name 
"Simoramic". 
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4. 



The PPM Spectral Model 



Pulse -position modulation (PPM) is a form of pulse time modu- 
lation in which the modulating process causes the position (in time) of 
the pulse to vary relative to its unmodulated time of occurrence. The 
relation of the pulse to the modulation process is normally determined 
by the sampling technique. The two more common types of sampling 
are "natural sampling and "uniform" sampling [2] . For natural sam- 
pling, the occurrence of the pulse is determined by the instant at which 
the leading edge of a sawtooth sampling wave and the modulating (mes- 
sage) voltage are equal. 




Figure 4-1. Generatipn of PPM by Natural Sampling. 

When uniform sampling is used, the modulating waveform is 
sampled at equally spaced time intervals. Each pulse in the PPM wave- 
form is then delayed in time from its unmodulated position by a time e 
which is proportional to the sampled voltage. 
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message waveform x(t) 




If one desires to relate a specific modulating process to the re- 
sultant power spectral density of the PPM waveform, it is necessary 
that the pulse positioning process e be expressable in terms of the mod“ 
ulating process x. This was done for the case of uniform sampling; 
however, the attempt to do this for the natural sampling case proved 
fruitless. The resulting expression for the natural sampling case was 
in the form of an unwieldy series of expansion which could not be reduced 
to a closed form nor was of any value in its series form. Therefore, 
the spectral model development will be based on a single channel PPM 
system using uniform sampling. 

The zero order hold is employed by the uniform sampling tech- 
nique just described. The model can be adapted for use with other sam- 
pling schemes. For a Gaussian process, one can use a higher order 
hold in the model, as the model will be developed to handle the case 

An example would be a technique which requires that the value 
of €, be linearly related to the average of the previous N samples of 
x(t). 
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for which € , e are correlated Gaussian bivariate random variables 

m n 

for a finite m-n. The correlation for this hypothetical case can be in- 
troduced by the sampling technique, or it may be directly related to the 
correlation of the sampled values of the modulating process, or both. 

Let the modulating waveform x(t) vary from zero to A volts. 
Define the maximum allowable deviation of the pulse from its unmodula- 
ted position as T /L seconds. The time shift e, is then 
o k 



T T 

e k = Ha x(kTo/ LA 



The pulse time shift e is thus linearly related to the sampled amplitude 

K 

of the message waveform x(kT^). The variance of is also linearly 
related to the variance of the message waveform samples x~ or 




T T 2 

-eI w x} • or 



2 

a 



€ 




The standard deviation (or RMS deviation from the mean) of the pulse 

time shift e, is 
k 



a 

e 



T 

— a 

LA x 



A more meaningful expression in determining what the actual 



RMS deviation is relative to a normalized interval will be defined as CJ . 

n 



a a 

CT = _£ = _JL 

n T LA 
o 
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A general expression for the power spectral density of a PPM 
waveform will be shown. The process x(t) producing the modulation 
will be considered stationary. The development of the equations will 
follow very closely that of Middleton [8] . 

A pulse function U(t - nT^ - e) shall be defined as having any 
desired shape (with maximum value of one) and existing only for a dura- 
tion u beginning at time t = nT^ ^ e n * Outside of this specified inter- 

i 

val, the pulse function is zero. Define y^~ as the peak amplitude of the 
th 

i — member of the ensemble of PPM waveforms in the time interval 

nT +e < t < nT + c + u. A member of the truncated ensemble 
on on 

y^.(t) may now be written: 

N 

*N*W “ I 0®“ " " T o - O 
n=~N 

The autocorrelation function of the PPM waveform 



i, , lim i. v 
y ^ " N-oo y N ^ 



is R (T). 



rV) = lim 1 r 



T/2 



T-oo T J y N (t)y N (t + T)dt 



-T/2 



lim 

N -*oo UN + 1 )T 



N 



hrr I 

° m, n=~N 



i l 
m n 



pOO 

U(t - mT 



-oo 



e l )U(t + T - n T “6 1 )dt 
o m on 



Now let X^(f) be the Fourier transform of U(t). 
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, CO 



X TT (f) = J U(t)e" jlUt dt 



U' 



-oo 



or 



puu 

J U(t -a)e" jWt dt = X TT (f)e" JUja 



U v 



Also 



U(t) = J X TT (f)e J ~ k df 



Jwt 



-oo 



U x 



or 



U(t - P) = J X u (f)e jUj(t " P) df 

~oo 



Letting a = nT +0 - T and P = mT + € in the above equations, 

on o m 

have 



rV) = 



liin 



N-oo (2N + 1)T 



• • P 00 

X u(t - 

n = -N ~OQ 



a) 



0 m , n=“N 



n 00 



-00 



X u (f)e :il,j(t “ P) dfdt 



lim _ 
N -00 (2N + 



1 V i if 00 ,, 

f 1 )T L y n y m J X U 

° m, n=-N “°° 



(f)e 



-jiuP 



p 00 
'■too 



U(t - a)e ja,t dtdf 



lim 1 



N -00 (2N + 1 )T 



1 v»J 



OO 



m, n=-N 



00 



X u (f)e 



-juuP 



x u( f)V“df 



we 
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rV) = f J 



oo 



o -oo 



|x u <f)| 6 In-co (2N + 1) 



N 



V i i n 

/ y y e 

Li n m 

m, n=-N 



. juu{e -e 1 ) juu(n-m)T 
- - m o 



’] 



df 



Now is it necessary to determine the covariance function from the above. 
Recall that the covariance function is defined as 



K(T) = E 



{rV)} 



Additionally, recall that it is necessary to perform the ensemble aver- 
age prior to taking the time average of the autocorrelation function. 
Thus 



n oo 



K(t) = — 



o -oo 



|x TJ U,|V^ 



1 



liitr 



N - > oo (ZN + 1) 



N 

{ z 

m , n=-N 



. juu(e 1 -e X ) juu(n-m)T 
ii n m 

y y e e 

m 'n 



°} 



df 



Now observe three important facts. The expected value of a finite ser- 
ies is being taken, such that y vanish for m , n > N. Since 

m, n jiii 

the process is stationary, the statistical relationship of any two random 
variables displaced the same distance kT in time is the same for all 
such pairs of random variables. Finally, there are 2N + 1 such like 
processes in the double series. Define k = m - n in the series, and the 
covariance function may be ^rew'ritt en 
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' 




K(t) = 



oo 



T j IV f) ' Vi " T 

O "OO 



lim 

N-oo 



N 

I 

k=-N 



ju>(e -e^ 



E ( y l y 2 e 2 1 } 



juukT 



df 



t -t =kT 
2 1 o 



Now that the mean of the ensemble has been taken, one may apply the 
limit, and 



p°° r i ? r 

K<T) = J T l x u (f) l l 

"OO O 1„ ^ 



: { y l' 



E iyiy 2 e 



k=-oo 

juj(e z -e 1 ) 



} 



juukT 



t_-t=kT 
2 1 o 



-jwT 
e df 



K(T) is now in the form of an inverse Fourier transform; thus the 
Wiener-Khintchine Theorem for two sided density functions may be uti- 
lized at this point. Let S l (f) be the two sided power spectral density. 

1 7 v? r jo)(e -e ) jwkT 

S '(D = f Ix^fll £ E{ yi y 2 e } e 

° k=-oo kT 



-oo < f < oo 



Later, the power spectral density shall be defined for positive fre- 
quencies only (one sided), and since S 1 ( f ) is an even function, the one 

sided power spectral density is S(f) = [2S , (f)]u (f), where u (f) is 

s s 

the unit step function. The unit step function is defined as 
0 , f < 0 



u (f) = 1/2 , f = 0 
s 

1 , f > 0 
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Define S^'(f) as the two sided spectral density for the typical sampling 



pulse, or 



lx (f)| 

V (f) = “ T * -co *f *oo . 



S 1 (f) is a continuous even function of frequency. The one sided spec- 
tral density for the typical sampling pulse would then be 



s u (f) = 2 S u ,(f) = 



2 iXytf) l' 



, f * 0 



= 0 



, f < o . 



Thus, the power spectral density may be rewritten 

juukT 



oo 



s '( f ) = Sy-if) £ E {yiy 2 e } 

k= -oo kT 



, -oo ^ f ^ oo 



The spectral model to be used assumes that the y^ are indepen- 
dent from sample to sample, and independent from the £ process. 

Thus, the equation that is used for the generalized spectral model is 



S^f) = S^f) 



oo 



vW I 

k = -oo 

k^O 



E|e } 

k W'l 



, -oo ^ f ^ oo . 



Observe that the model can also be used for pulse -amplitude modulation 

th 

(PAM) in which the y^ represent the voltage level of the k — sample of 
the message and the pulse is subjected to a random jitter (in time), 
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possibly due to equipment temperature fluctuations or varying transmis- 
sion path length. The y^_ would of necessity be independent random 

variables for the generalized spectral model just developed. The e , 

.K 

however, as will be shown in the further development of this PPM 
model, can either be independent, or correlated random variables with 
a Gaussian distribution. 
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1 



I 

! 






5 . 



A Discussion of the Statistics in the Spectral Model 
As seen in the development of the spectral model, certain sta- 
tistics of the two processes e and y, together with the pulse shape, in 
essence determine the power spectral density,, Since independence was 
assumed, the statistics of the pulse amplitude y are quite simply de- 

K 

termined, as all that is involved is the mean of the square of the pro- 

~ _2 

cess (y ) and the square of the mean (y ). The statistics of the e pro- 
cess are not as simple, especially since independence of e will not be 
assumed for all cases. The statistics of the e process are directly re- 
lated to the message process x(t). It has already been shown that the 
variance and the standard deviation of the two processes are linearly 
related. Furthermore, the probability density functions (pdf) of the two 
processes are of the same form if one considers continuous distribu- 
tions only (i. e. , if x is Gaussian, e is also Gaussian). Specifically , 

if p (x) is the pdf for x^ then 

x k 




To show that this is true, recall that e = (0 /a )x . Now 

e x 



P[ e < e] = P 
k 



a e 

X 



X. 



k a 



But 



a e 

. X 

X < 

k a 



p a 

= i 

-oo 



a e/a 
x e 



p (x)dx 

x k 



Since the pdf is the derivative of the distribution function, then 
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P- (£) 

£ k 



_d_ 

de 



per e/a 

X 0 

J 



-CO 



p (x)dx 
X k 



from which one can obtain the desired result 

a /a e N 

k G k \ e , 

To avoid carrying unnecessary terms in equations already complex with 
terminology, a time axis shift of the PPM waveform will be assumed 
such that the mean of the process 0 is zero. That this simplification 
may be performed with no effect on the resultant power spectra is in- 
tuitively obvious, but nevertheless has been previously verified by the 
writer in each case considered. 

To evaluate the statistics of the process 0 that appear in the 
equation for the power spectral density, one must determine the joint 
characteristic function G(j(JUj, juu^) °f a bivariate process. The joint 
characteristic function is defined as 

j ( Vi + V 2 >H r MVhh 

G(jO) j , jo> 2 ) = E|e | = E|e j . 



This joint characteristic function will be evaluated for the Gaussian 
bivariate case. For the remainder of the cases included in the compu- 
ter model, statistical independence from interval to interval of the 0 
process was assumed. For the independent case, the problem reduces 
to the determination of the square of the absolute value of the charac- 
teristic function of 0. 
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{< 



E e 



juu(e -ej). 



}:= |E {e JUJC } | 2 = |G £ (jaj) | Z 

The equation for the power spectral density now reduces to 

’ — oo 

y 2 + y 2 |G £ (oj)| 2 £ e ja,kT < 



S' (f) =S u '(f) 



k = -oo 
ky^O 



- oo ^ f ^ oo 



But 



oo oo 

>kT c J_ 

T 

-oo o -oo 



I 



H-f) 

-oo \ O / 



Recall that S (f) = 2S'(f) for f ^ 0 and S(f) = 2S’(f)u (f). Then 
U • s 



U v 



— 


r °° / \i 




2 _2 . , , , 2 


, 1 Vi/, k \ 




y +y l G (w)l 


- , + ¥- L^V-t) 


\ 


l 


0 

II 

0 

/ 

0 


I 



u (f) 

s 



for the case in which £ and y are independent from sampling interval 

KL .K. 

to sampling interval. As mentioned previously, only positive frequen- 
cies are considered. 
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5A. The Message as a Gaussian Process 

The purpose of this subsection (and the following subsections) is 
to fit a specific process into the generalized spectral model. The pro- 
cess considered may either be the message process x or the € process, 
as the similarity of the two processes has been previously shown. For 
the Gaussian case, the development will allow for a statistical correla- 
tion to exist between I adjacent samples of the process. More speci- 
fically, this development will allow for the spectral model to handle 
Gaussian processes for which the autocorrelation function R(T) exists 

for x(t) up to I T | = IT , and R(T) = 0 for I T I > IT „ This is practi- 

o o 

cal because a purely random process yields an autocorrelation function 
that approaches zero as T increases [l]. 

The expected value 



E {e >(C Z- €l) } 



kT 



is the joint characteristic function of e^, £ ^ , for a given kT = t^ - t^, 
The joint characteristic function G( j(JU^ , juu^ ) as previously defined fits 
the case now under development when (JU = -(JU . It was decided to 

-L Ld 

carryout the general development of G( juu ^ , jw^), substituting into the 
final result w = (JU ^ = *00 j. Thus, 



jw(e -e ) 



{ e 2 1 } ~ E {‘ 



j(w e +w e ) 

, 11 z 2 e u 



} = G(ja) 1 ,juj 2 ) 



The expected value is determined as 
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->oo poo (j(u 1 e 1 +juj 2 e 2 ) 



pUj, e 2 )de 1 de 2 



-oo -oo 



The bivariate Gaussian pdf is 

1 



p(e 1> e 2 ) = 



[ - 1 / 2ct 2 ( 1 _ p 2 ) ] (e 1 2 -Zpe 1 e 2 +e 2 2 ) 



2TTG 2, yji~- p 2 



where 



\Jl = 0 and O = (J = O 
e e i e 2 



For simplicity, substitute 



Sj = X 



e 2 = * 



20 = a 



2 



p = b 



Then 



G(jou 1 juu 2 ) 



Trab 



2 2 

oo j(D / ab ) 

e 



-oo 



oo 

r 



-oo 



2 W , 2 






-(1/ab )(X -2f*Xi|r) J 1 



dXdijr 



But 



p oo 
-oo 



-(l/ab 2 )(X 2 -2pX^) j V 
i e dX 



2.2, 2 o oo 

= e P * /ab , e 



. p UU 

J-oo 



-(l/ab 2 )[(X-pt) 2 ] 

e dA. 



which is in the form of a constant multiplying the characteristic function 

2 2 

of a single Gaussian process with mean | d = p\|f and a = ab / 2. 



The characteristic function of a single Gaussian process may easily be 
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determined (with the aid of tables of integrals), and the last integral is 



simply 



b Var’’ 



2 . 2 . ,2 
p i|r /ab 



2 2 

jpilroUj -(ab oUj / 4) 



The joint characteristic function then becomes 

(p 2 -l)/(ab 2 )* 2 



rr . , 1 -{ab 2 u) 2 /4) o0 ° 

G(j “V a ’2 ) = vST 6 1 



Define 



I 



dt 



oo 



UU 3 = UJ 2 + p'JUj 
2 1 ( ab 



2 V 2 

P - 1 / 

and this integral also takes the form of the single characteristic function 
of the process i|; , with the mean u = 0 



G(jo ll ,j«, 2 ) = 0^ 



rr 


l p°° -(!/2a^ 2 )t 2 jaJ 3 4f 


k e 


/ — e e dll; 

L a t ^ -oo J 



• 2 -(ab 2 uu 1 2 /4) 

= a Wr e 



•<° t V /2) 



By re substituting the values for a^, a, b, and one obtains 

‘Z + ' U 2 



z z z 

■(o / 2)(uj 1 + 2puj uu +tju ) 

G 1 1 Z Z 



G(ja)^ 2 ) = e 

Recall that for this development, oj = (JD^ = Thus 3 



G(j<ju 1 jt« 2 ) = G(ou) = e 



2 2 .- , 
-cr uj (1-p) 



Applying the above to the spectral model for a given value of kT^, we 
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have 



E f jtu(e 2 -e l 

J kT =t -t. 
o 2 1 



-a V[l- P (kT )] 

G O 



That the expression is a function of the autocorrelation of the process e 
(or x) may be determined by observing that p is the correlation coeffi- 
cient for a given kT in the development. Considering all kT^, p(kT^) 
is the normalized autocorrelation function. 



p(T) = 



R(T) 

R(0) 



Now let R(T) = 0 for k > I. The power spectral density for this 



case is then 



( 



r i 



S'(f) = S^f^y 2 + y z 



2 2 



I 



C uu [l-p(k)] juukT 



k=-I 

k^O 

-a V/- ( i, + 1) j»kT ” jiDkT 

+ e ( 2, e + 2, ' 



k= -oo 



k=(I+l) 



But p(kT ) = p(-kT ), as R(t) is an even function, and 
o o 



S'(f) = Sy'tf) ( y 2 + y 2 



* / -a e 2 uu 2 [l- P (k)] -cy* 2 ’ 

2 1 l e 

k=l V 



— e 



cos uukT 



oo 



juukT 



2 2, 

-o’ w / 

te £ (- 1+ k L e ) 



and the spectral model for the Gaussian case (for f ^ 0) becomes 



S(f) = S^f) { y 2 + y 2 



T / 2 2 r . „ 2 2 k 

i, / -O 11) [ 1 - p(k) j -CT 

^ 1 G G 



I 



- e 



cos uukT 



L k=l 



4- e 



2 2 

-a uu 

G 



oo 



o 1 ,_ 



k = 0 



u (f) 
s 
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5B. The Message as a Uniform Process 

The specific spectral model will be developed for the case in 
which the message determining the modulation of the PPM waveform is 
a uniform process. Statistical independence of both the e and the y 
processes will be assumed. For this assumption, the spectral density 
was shown to be 



S(f) = S u (f){y 2 + y 2 |G e (uu)| 2 



oo 



1 + 



f I 

o k = 0 



6 f 



k (f) 



The problem is therefore reduced to that of determining the character- 



istic function G^(uu) of a uniform process. 



G (uu) = e UJG p (e)de 

£ d ^-i 

-oo k 



As shown previously, if x is uniformly distributed, then c is also uni- 
formly distributed. For the case of uniform distribution of x, the total 

possible deviation e(max) allowed is T /L seconds (or ±(T /2L») from 

o o 

the zero mean assumed). Therefore, 



T 



P e, (e) T ’ ~ 2L ^ € ^ 2L, 

k o 



Thus , 



G (uu) 
c 



T /2L T 

o L juue 

T~ e 

-(r / 2Li) o 

o 



de 



o r 



uuT 



sm 



G («j) 
e 



2L 



uuT 



2L 
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The power spectral density for the uniform case then becomes 



S(f) 



S l / f) \ y2 + Y 2 





o 
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5C. The Message as an Exponential Process 

The specific spectral model will be developed for the case in 
which the message determining the modulation of the PPM waveform is 
an exponential process. Of interest for this case is that the speech 
process is approximately exponentially distributed, especially in the 
region of medium to large amplitudes [1] . 

For the exponential case, statistical independence will again be 
assumed for both the 0 and y processes. As was the case then in sec- 
tion 5B, the problem is reduced to determining the characteristic 
function of the process 0 . 

The exponential probability density function is 



p (e) 
e k 



1 : (V2/a e )|e| 

v^T e 



from which one can determine 



|J = E(e) = 0 

0 

Var(s) = E(e 2 ) - H 2 = E(e 2 ) = a 2 

0 0 



Define 




0 



Then 



P (e) 

e k 



a -a e 

— e 1 
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The characteristic function is 



pOO . 

G (uu) = e JUj€ p (e)dc 
J -oo e k 



a 

2 



0 e €(j«,+ a ) d£ + ( ,C ° e eO»-a) de 

0 



r 

-oo 



or 

C e ( *> = t£t 

au + a 

The power spectral density of the PPM waveform for the exponential 
case is then 



S( f > - v f V 2 + y 2 , ■ i ? 

(id + a ) 



- 1 + 



f !•('-£)] 

O ' O ' -I 



U (f) 
s 
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6 . 



The Computer Model 



The PPM spectral model was developed in section four, and the 
model was adapted to specific modulation processes in section five. 

The resultant equations are quite formidable if manual calculation is re- 
quired to determine the power spectral density for various cases. For- 
tunately, the tedium of manual calculations of such proportions is not 
necessary nowadays, as the digital computer handles such tasks with 
great facility. A computer model will now be developed that will incor- 
po rate all cases previously discussed, plus the unmodulated case. 

The computer model developed relates the pulse positioning pro- 
cess 0 in terms of the originalmaodulating process. This allows the 
user to have a more direct insight as to the effects of the modulating 
process on the power spectral density of the PPM waveform. Most of 
the parameters of the model are variable at the discretion of the user. 

An exception is the number of points of continuous spectra calculated 
between each possible discrete frequency spectrum. Since five points 
in each such interval are adequate for a graphic plot, this quantity has 
been fixed at five to avoid possible confusion by future users of the 
model. 

It should be noted here that the following development of the com- 
puter model may not at times appear to be the simplest approach. The 
approach taken generally yields simplification of computer Fortran oper- 
ations to be performed, such as eliminating the calling of a library 
function for performing an exponentation when a simple multiplying 
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operation will suffice. It is admitted, however, that a few of the "sim- 
plifying 11 techniques employed did not necessarily make the program 
more efficient. The development will include these also, as they have 
become a part of the final computer model. 

While no extreme efforts were applied toward making the com- 
puter program efficient, this aspect was certainly kept in mind while 
the program was being developed. 

The program was written in Fortran IV language for the IBM 
7 0 40 at Sylvania 1 s Electronic Defense Laboratories (EDL). With some 
obvious card changes (for each library function used), the program can 
be adapted to the Fortran language of the CDC 1604 at the Naval Post- 
graduate School. 

The total run time (including card reading and program printout) 
for computing four separate cases in one program run from a binary 
deck was timed at two minutes on the 7040. 

From the previous developments, it was decided to base the 
computer model on the mathematical model expressed in the following 
form for frequencies greater than zero. The case for d-c (zero fre- 
quency) will be developed separately. 

S ( f ) “ S u (f){y 2 + y 2 [Sum(f) - | 0^(00) | 2 J j- 

oo ^ 

+ s u( f >{7 2 l G € (»>l 2 f } 

° k=l ° 

The first half of the expression is the continuous spectrum and the 
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second half is the spectrum concentrated at discrete frequencies. The 



term Sum(f) is defined as 

I /„ .,,.0 ... v 2 

Sum(f) i 2 ^ 
k=l 




The term Sum(f) is included in the model so that the Gaussian correla- 
ted process may be handled. Note that when Sum(f) equals zero, the 
model reduces to that derived for the case in which all processes were 

assumed independent. The term G ((ju) is simply the characteristic 

S 

function of the independent random process e, and fits as shown into 
the chosen representation even when the Gaussian correlated case is 
considered. That this is true may be determined from the fact that the 
characteristic function for the independent Gaussian case (mean zero) 
is the square root of that determined for the Gaussian bivariate case 
when rho (p) is zero, or 

“[(<* U>) 2 /?] 

G (uo) = e 
0 



(for Gaus sian proces s with p = 0). 

All terms in the model have now been adequately defined except 
for the power spectral density of the "typical sampling pulse 11 . Recall 
in a previous development that this pulse could be of any shape with a 
maximum amplitude of one. Also recall that this amplitude of one does 
not restrict the amplitude of the pulses in the PPM waveform, as the 
PPM pulse amplitudes are determined by the y process. Thus, the 
"typical sampling pulse 11 is simply a description of the shape of the 
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pulses of the PPM pulse train. For this computer model, a rectangular 
pulse will be assumed, the pulse being of duration u and of unit ampli- 
tude. 



bH 



Figure 6~1. The Typical Sampling Pulse for the Computer Model. 



As previously shown, 



s u (f) = ^ |x TT (f )| 2 . 



IF 



Here 



oo 

X u (f) = I U(t)e“ JUJt dt = 

-CD 



p u/2 
“(u/2) 



(lle-j^dt. 



from which one may determine 



- 2 / sin — 

S (f) = — 2 

IT 1 T V urn 
o 



To simplify the computer program (so that some of the time consuming 
operations are minimized), the above equation may be rewritten as 



c if\ (1 cos (juu) 

S U (f) T .2 

o (am) 



The discrete frequencies at which S(f) is to be computed are 



. n 2 irn . . „ „ . . 

f = , or oi = y-qr - , 0 ^ n ^ 2 00 



n 5T 



n 5T 
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To caiculate S(0) 5 one must apply the value of the unit step 
function (one -half) at zero frequency (the d-c term). Substituting 03 = 0 
into the equations developed thus far 5 one obtains for all cases 



S(0) 




where 



V 0) 

z 



2 

u 

T 



o 

The sequence of computation of the power spectral density may 
be determined from the flow chart (Figure 6-4). 

The computer model will determine two hundred values of S(f) 
up to and including f = 40/T^. This includes all 40 possible values of 
S(f) concentrated at discrete frequencies (delta functions) within the fre- 
quency range considered. 

A list of the significant symbols used in the program and their 
relationship to the terms in the mathematical model is shown in Table 
6- 1. The computer program itself with one page of output is included 
as Figures 6-2 and 6-3. 
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Table o-l. 



Computer Program Symbol Identification. 



Symbol 


Equivalent Symbol In Mathematical Development or 
Symbol Definition 


A 


A 


ALFA2 


CL 2 (used for exponential case) 


CFSQ(N) 


|G e (UJ;| Z 


CONT(N) 


Value of continuous spectrum at f = N/5T 

o 


DELFN 


Value of discrete spectrum component 


FL 


Li 


F(N) 


Frequency; F(N) = N/5T^ cps 


I 


1 (upper summation index for Gaussian case) 


J 


None; this index specifies modulation process 8 See 
comments in Figure 6~2« 


NR 


Number of separate cases to be computed this program 
run 


RHO{K) 


p(k) 


SIGMA 


a 

X 


S(N) 


S(f); Power spectral density 


SU(N) 


S u (fi 


SUM(N) 


Summation for Gaussian correlated case 


T 


T 

o 


TEMP(K) 


Partial value of summation of SUM(N) 


U 


U 


YBSQ 


2 

y 


YSQB 


~z 

y 
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PROGRAM PPM 

THIS PROGRAM CALCULATES THE POWER SPECTRAL OENSITY OF A PPM 
WAVEFORM THAT HAS AS ITS MODULATING SOURCE ONE OF THE FOLLOWING 
EITHER AN INDEPENDENT GAUSSIAN PROCESS OR A GAUSSIAN PROCESS 
WITH UP TO TEN CORRELATION VALUESIFOR TAU=T , 2T , . .. , 10T ) , AN 

“Independent uniform process, or an independent exponential process 

J=1 FOR GAUSSIAN, *2 FOR UNI FORM, =3 ( UNMODULATED ) »=4 FOR EXPONENTIAL 
FORMAT 8 STATEMENT ALLOWS FOR READING RECORD IN FOLLOWING ORDER 
U»T,YBSQ»YSQB,SIGMA,A,FL,J»I» AND 10 VALUES OF RHO 
■DIMENSION' RHO (TO) , SUM (200), CFSQI 200! , SU< 200), S(200) .TEMPI 10) , 

1FI200), CONTI200) 

“ 6 FORMAT (4H F= F10:3',“6H SIFI* E12.4, 9H CONT(F)^~EI2.41 

7 FORMAT I/5X, l HU, 11X,1HT,9X, 4HY8SQ.8X, 4HYSQB,6X,5HSIGMA,6X, 1HA.10X, 
I 2HFL , 4X , 29H 2ND ROW IS J,T» RHOI K= l , TO) /) 

8 FORMAT (7F11. 6/(213, 10F6. 3) ) 

9 FORMAT (AH F* F10.3, 6H S ( F ) « E12. 4) 

PI= 3.16159265 
NR - 4 

300 PRINT 7 

REAO 8,0 , T , Y B S Q ,YSQB, SIG M A , A , FL , J, I ," ( R HO ( K ) , K = 1 , 10 ) 

PRINT 8, U,T, YBSQ.YSQB, SIGMA, A,FL, J, I, (RHO(K) ,K»1, 10) 

00 301 N«l","2tf0" 

301 SUM (N ) = 0.0 

GO TO ( 1,2, 3,4) ,J 

l C= (0.4*PI*SIGMA/(FL*A) )**2 
DO TlN=l ,200 
FL02 = N 
CC=C» FL 02* F L02 
IF (CC.LT.23.0) GO TO 12 
CFSQ(N) = 0,0 
GO TO 11 

12 CFSQ(N) * 1.0/ ( EXP (00) )'' 

11 CONTINUE 

IF Tf.EQ. 0) GO TO 100 
DO 18 N=l , 200 

FLCT3 U 

TEM - 0.0 
DO 15 K=1 , 1 
FL04 = K 

CD =* C4FL03*FL03*( l.O-RHO(K) ) 

IF (CO. LT. 23.0) GO TO 14 
TEMPT K) * OiU 
GO TO 15 

14 TEMPI K) = (1.0/(EXP(CD) )~CFSQ(N) )*( COS (0.4*PI*FL03*FL04) ) 

15 TEM - TEM*TEMP ( K ) 

18 SUM (NT * 2.0*TEM 

GO TO 100 
"Z PU * 0.Z4FT7Ft 
DO 25 N=l , 200 
FL06-N 
PUN=PU*FL06 

25 CFSQ( N)=( ( SIN(PUN) )/PUN)** 2 
GO TO 100 
3 DO 35 N=l, 200 
35 CFSQ(N) * 1.0 



Figure 6-2. Source Program for PROGRAM PPM, Page 1 
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GO TO 100 

4 ALFA2 = 2.0MFl*A/(T*SIGMA) )**2 
PE = ( 0 . 4* P I / T ) * * 2 
DO 45 N= 1 » 200 
FL05 = N 

45 CFS0( N)=(ALFA2/(PE*FLG5*FL05 ♦ ALF A2 ) ) * * 2 
100 AA = 4.0*U*U/T 
B= 0.4*PI*U/T 
DF= 0.2/T 
DO 110 N= 1 , 200 
FL01 =N 

F ( N ) = FLO 1 * DF 
G = B *FL0 1 

110 SU(N) = AA*( l.O-COS(G) )/(G*G) 

SUZERO = AA/4.0 

SZERO = SUZERO* ( YSQB -YBSQ* ( 1.0-1 .0/T ) ) 
FZERO= 0.0 
DO 120 N = 1 t 200 

120 S (N ) = SU(N)*(YSQB+YBSQ*(SUM(N)-CFSQ(N)>) 
DO 125 N=5 f 200 * 5 
CONT(N) = SIN) 

DELFN = SU(N)*YBSQ*CFSQ(N) *1.0/T 
125 SIN) = CONT(N)+DELFN 
PRINT 9, FZERO, SZERO 
DO 210 K~- 1,196,5 
K3 = K+ 3 

~ DO" 200 N=k,K3 
200 PRINT 9, F ( N ) , S ( N ) 

N=kV4 

210 PRINT 6 , F ( N ) »S(N) ,CONT(N) 

NR = NR-' 1 

IF (NR.GT.O) GO TO 300 

C" THIS RUN FOR 10 KC SAMPLING. U = 0.1T 
STOP 
END 



Figure 6-2, Page 2 
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Figure 6-3. A Sample Output from PROGRAM PPM 
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7. Computed Power Spectral Densities, Results 

This section contains the displays of the power spectral densities 
computed for certain combinations of pulse widths, modulation processes 
and varying degrees of pulse deviation (or RMS deviation from the mean 
position). The legend accompanying each graph indicates the pertinent 
data regarding the process for which that power spectral density was 
computed. 

The first graph (Figure 7-4) is the power spectrum of an unmo- 
dulated pulse train. Observe that with no modulation of the pulse train, 
there is no continuous component in the power spectrum. 

Three different normalized autocorrelation functions were used 
for the Gaussian case. For the process which has very little correlation 
between samples, the approximation used was the linear autocorrelation 
function indicated by Figure 7-1. 




Figure 7-1. Limited Cor relation Between Samples. 

An approximation to the autocorrelation function of a process with a high 
degree of correlation existing for several adjacent samples is shown in 
Figure 7-2. 
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p{T) 




Figure. 7-2. Samples Highly Co r related. 



An approximation to the autocorrelation existing for band-limited Gaussian 
noise (0~W cps) is shown in Figure 7-3. In Figure 7-3, the Gaussian 
noise is assumed to have been sampled at ten times the highest frequency 
present^ or at 1 OW samples per second* 




Figure 7-3. Approximate Autocorrelation for Oversampled 
Band-limited Gaussian Noise* 

Two examples are shown for which the pulse amplitudes are not 
held constant (Rayleigh distributed), in addition to being pulse position 
modulated* It may be observed that in these cases the continuous spec- 
trum begins with an appreciable value at zero frequency, which is a 
characteristic of pulse amplitude modulated waveforms. 
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To allow for an adequate display of the continuous spectrum, a 
scale facto ? was chosen that does not permit the showing of the discrete 
power spectrum to scale* The values of the discrete power spectral 
lines are numerically indicated for each line that exceeds the scale of 
the graph* 
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-10, Spectrum for Exponential Case #3. 
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-16. Spectrum for Gaussian Case #6 
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Spec ral Analysis of Actual PPM Waveforms 



A comparison of the results from the spectral model was made 
with the results determined experimentally. The objective of this aspect 
of the investigation was to determine if short time interval analysis 
would yield a reasonable estimate of the power spectral density of PPM 
wave fo rms. 

In order to generate PPM waveforms comparable to those asso- 
ciated with the spectral models a pulse position modulator was designed 
and built* The circuit and performance characteristics of the modula- 
tor are included in Appendix A c The unmodulated pulse train generated 
has a pulse repetition rate of one kilocycle s and a pulse width of 195 
micro seconds. 

The spectrum analyzer (periodogram calculator) used for the 
experimental measurements is known by the trade name M Spectran'\ 

The Spectran is one of the parallel filter bank types of analyzers des- 
cribed in section 3. It contains 480 parallel filters s and has a resolution 
of 25 cycles per second (cps)* The analysis time is 18 milliseconds. 

The first two frequency bands of the Spectran were used for the PPM 
analysis, covering 20 to 3852 and 3852 to 7684 cps s respectively. The 
use of two bands accounts for what appears to be a slight discontinuity in 
the pictures to be presented later in this section. 

The Spectran was calibrated by inserting a sine wave input of an 
arbitrary frequency (two kilocycles) and recording the analyzer response 
for measured power inputs* The input power was considered the square 
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of the RMS voltage input* which is equivalent to the spectral model defi- 
nition of power Into a one ohm resistor. The resulting calibration curve 
is included as Figure 8-1* Unfortunately, the calibration curve is only 
accurate for spectra at discrete frequencies (delta functions) as a conse- 
quence of having been based upon data taken for discrete frequency in- 
puts. The calibration curve will not yield an accurate value of the con- 
tinuous spectrum at a specific frequency f^. Each of the 480 '‘spectral 
windows" of the Spectran has a bandwidth of about 10 cps. Therefore, 
when a continuous spectrum is measured, the measurement obtained at 
a frequency f^ is actually the proportionate power of all frequencies that 
are passed through the spectral window that passes the frequency f^. 
Further^ the noise level of the measuring equipment prohibited accurate 
measurements of the continuous portion of the spectrum. Consequently, 
for the experimental results, a determination of the exact values of the 
continuous spectrum must be attempted by a different procedure. The 
general shape of the continuous spectrum was determined by the Spectran 
analysis, however and is displayed in Figure 8-3. A discussion of the 
background related to the photographs in Figure 8-3 is considered neces- 
sary in order that they may be properly interpreted* As predicted by 
the computer model, the value of the power spectral density at any fre- 
quency of the continuous spectrum is several orders of magnitude less 
than that of the larger discrete components. In addition, the predicted 
average levels of the continuous spectrum is about the same as the noise 
level output of the measuring equipment (the Spectran). Thus, if one 
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visually subtracts the noise level and the effects of the discrete compo 
nent as determined from the top picture in Figure 8-3, then a fair esti. 
mate of the shape of the continuous spectrum may be gleaned from the 
bottom picture. 
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Figure 8-1. Spectran Calibration Curve 
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Figure 8-2. An Experimental Estimate of the Discrete Power Spectral Density. 






Spectrum of the 
unmodulated PPM 
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Figure 8-3. An Experimental Estimate of the Continuous Power Spectral Density. 




In order to measure the value of the continuous power spectrum 

at d screte frequencies, a wave analyzer with a tunable 6 cycle bandwidth 

was employed* The PPM waveform, modulated by a random noise 

(assumed Gaussian), was used as the input to the wave analyzer. True 

RMS voltage measurements of the analyzer response were made at 100 

2 

cycle increments,, The block diagram of the experimental setup is 
shown in Figure 8-4* 




Figure 8~4. Technique for Experimental Measurement of 
Continuous Spectrum. 



Assuming an ideal rectangular response of unit amplitude for the 
filter, the total power to the RMS voltmeter is 

nf +3 

Power = I ° S(f)df , 

J f “3 
o 

where is the center frequency of the filter. It shall be assumed that 

The analyzer is the HP-302A, which has an equivalent noise 
bandwidth of 6 cycles with the center frequency adjustable throughout the 
audio frequency range* 

Z 

True RMS voltage measurements were made with a "Hi-damped” 
Fluke Model 910 - A voltmeter, which is a power sensitive instrument 
calibrated in volts RMS. 
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there are no extreme changes in the power spectral density in any 6 cycle 
interval then S(f) may be approximated by a constant S(f Q ) e The meas- 
ured power Is then approximately 

Measured Power ~ 6S(f ) . 

o 



S(f ) may now be determined from the square of the true RMS voltage 
measurement, which s the measured power (to a one ohm resistor). 



V 

s(f ) = — 

o 



2 

RMS 
~ 6 



A plot of the power spectral density measured in the manner just 

described is shown in Figure 8~5. The few inconsistent measurements 

recorded are attributed to this writer 5 s inability to ’’time average” a 

3 

slightly fluctuating voltmeter needle by visual observations. For com- 
parison purposes, the theoretical plot of power spectral density is 
included as the dashed line in Figure 8“ 5. 

The measurement of the discrete spectrum was a simpler matter, 
as output levels well above the equipment noise levels were being meas- 
ured, As previously mentioned, the calibration chart applies to the 
measured results of the discrete spectra, from which may be determined 
a reasonably accurate measured power level. Unfortunately, the smooth- 
ing of the Spectran "readout” is not perfect, resulting in the occurrence 



For a discussion of the variance of the output of a square law 
device (the same type as the measurement fluctuations encountered), see 
page 256, Davenport, W. B. and W, L. Root, An Introduction to the 
Theory of Random Signals and Noise, McGraw Hill, 19 58. 
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of peaks and valleys when one sweeps across the spectrum with a con- 
stant level nput s'gnal. Thus 5 even the discrete frequency power meas- 
urements w 11 have some equipment measuring error introduced. 

The discrete spectrum displays are shown in Figure 8-2. From 
F gure 8 2 and theoretical computations ? an experimental versus theore° 
tical comparison was made and is included in Table 8-1. The PPM de- 
vice is that described in Appendix A modulated with a wide band Gaussian 
noise source. 
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Table 3-1. Experimental and Theoretical Values of 
D screte Power Spectrum. 
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Discuss on of the Results and Conclusions 

A reasonable question to be answered now is "What specific sta- 
tistics of the modulation process contribute to the shape of the power 
spectral density ?" This is one question that is answered in this section. 
Other items discussed in this section are: 

1 0 The "general shape" of the power spectral density of PPM 
waveforms. 

2. The power spectral density as a function of the modulating 
process. 

3* The parameters that determine the relative amount of the 
total power that is contained in the discrete and continuous 
spectra. 

4. The effect of correlation between pulses in the pulse train 
with respect to their relative positions in time. 

5. Some general comments on the experimental determination 
of power spectral density. 

The shape of the power spectral densities of PPM waveforms 

2 

appears to be generally of the form sin x/x) if one considers only the 
frequencies beyond the first lobe. The portion beyond the first lobe has 
been observed to be mostly continuous spectrum for the modulation pro- 
cesses considered. The reason for this can be determined from the equa- 
tion upon which the computer model is based. The characteristic function 

1 

It will be assumed that all pulse waveforms discussed in this 
section consist of rectangular pulses except when otherwise stated. 
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for all modmat on processes considered "damps out" toward zero value 
as the frequency s increased, resulting in a power spectral density at 
higher frequencies that s essentially the same as the spectrum S^(f) of 
the typical sampling pul se* 

The frequency at which the first lobe is terminated (where the 
envelope of power spectral density first becomes zero) is the reciprocal 
of the pulse width* The shape of the first lobe appears to be quite depen- 
dent upon the modulating process* Figures 7-5, 7-8, and 7-14 are 
representative of the spectrum of PPM waveforms that differ in the modu- 
lating process i. e. , a uniform process, an exponential process, and a 
Gaussian process, respectively. The biggest difference in these power 
spectral densities is in the distribution of power in the discrete frequency 
components. ^ A close look at the Figures 7-5, 7-8, and 7-14 addition- 
ally reveals a distinct difference in the continuous spectrum within the 
first lobe* As was the case with the discrete spectrum, this difference 
may be attributed to the type of statistical process used to modulate the 
PPM waveform. 

The question asked at the beginning of this section will now be 
answered. For the cases in which statistical independence between 
samples is assumed the square of the absolute value of the character- 
istic function °s the only statistic of the modulating process that influ- 
ences the power spectral density of the PPM waveform. For all cases 

^ The frequencies at which the discrete components exist are 
multiples of the reciprocal of the sampling rate, or the average pulse 
repetition rate. 
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considered :,n this investigation, the absolute value of the characteristic 
function is s mply the Fourier transform of the probability density 
function of a process with zero mean, The process with a non-zero 
mean yields a characteristic function that is equal in magnitude to the 
characteristic function of the process with a zero mean, but with a phase 
difference. Taking the absolute value removes this phase difference, 
thus yielding the same results as for the process with zero mean. 

For the case in which statistical independence between samples 
is not assumed, the power spectral density is influenced by the joint 
characteristic function of the modulating process. The effects of corre- 
lated samples will be discussed in greater detail in the latter part of this 
section. 

The relative amount of the total power in the continuous spectrum 
of a PPM waveform with pulses of fixed duration and amplitude appears 
to be determined by the degree of fluctuation of the pulse shift. This ob- 
servation may be confirmed by comparing figures based upon like modu- 
lating processes but with different standard deviations (or different total 
pulse shift permitted in the uniformly distributed case). Comparing 
Figures 7-11 and 7-14, for which O' is 0. 05 and 0. 1, respectively, one 
may observe that an additional 21. 6 percent of the total power is included 
in the continuous spectrum for the greater fluctuating pulse shift. For 
the uniform case, the PPM fluctuations may be increased to the maximum 
level by permitting the maximum pulse excursion (delay) to be equal to 
the full sampling interval T^ less the pulse width u. This corresponds 
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u) in the PPM models developed. When this is done 



to L = T /(T 
o o 

the discrete spectrum becomes quite small and eventually disappears 
for very narrow pulse widths. The latter observation may be confirmed 
by direct calculations with the mathematical model; the former observa- 
tion may be confirmed from Figure 7~5. In Figure 7“ 5, the maximum 
pulse excursion within a sampling interval is permitted. The pulse width 
is 0. IT , and the aggregate power in the discrete components is reduced 
to less than one percent of the total power (excluding the d-c component). 

Another factor which influences the relative amount of total power 
contained in the continuous spectrum is the pulse width. For a given 
modulation process, as the pulse width is decreased, the percentage of 
the total power contained in the continuous spectrum increases. In con- 
trast to the considerations of the preceding paragraph, in which the total 
power of the PPM waveform is constant (independent of the fluctuations), 
if the pulse width is decreased the total power decreases unless a corre- 
sponding increase in pulse amplitude is made. A plot of the percent of 
total power in each component as a function of the pulse width is shown 
in Figure 9 ~ 1 for a single modulation process (the independent Gaussian 

case with CT - 0. 05). Data for this plot was obtained from the computer 
n 

model. 

The effect of statistical correlation between sample values of the 
modulation process tends to center the continuous spectrum about the 
discrete components. It is interesting to note that the relative amount of 
the total power in the discrete components remains essentially the same 
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for the correlated and the uncorrelated cases considered, yet the shape 
of the continuous spectrum is radically changed (see Figures 7-11, 12, 
and 13, for example)* However, in the limit as the correlation function 
of the modulation process approaches unity, the PPM waveform is essen 
tially modulated by a d~c level, and the power spectral density is all at 
discrete frequencies. In fact, the spectrum becomes that of an unmodu- 
lated pulse train (Figure 7~4) 0 

The correlated case is a practical consideration for typical PPM 
pulse trains encountered. For most PPM systems, there is time divi- 
sion multiplexing of several channels of data pulses, which are further 
interleaved with synchronizing (sync) pulses. The sync pulses are trans 
mitted at a periodic rate. The exact periodicity of the sync pulses re- 
sults in a contribution to the discrete components of the power spectral 
density* Further, the data pulses are normally statistically related to 
the sync pulses* The mean position of the data pulse is related to the 
sync pulse by some fixed interval for each process. The degree of cor- 
relation between data and sync pulses is set by the data pulse variance 
about its mean displacement* Finally, adjacent data pulses may be 
correlated, as was the case for the Gaussian correlated process of the 
simpler single channel PPM model used for this investigation. 

The verification of the theoretical results by taking experimental 
measurements was discussed in section 8, The measurement of the 
continuous spectrum of the power spectral density can be quite arduous. 
The difficulty is partly due to the fluctuations of the modulating process 
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(or the degree of modulation) varying in a probabilistic manner during 
the measuring interval* This effect, of course, can be minimized by 
averaging the measurements over a longer measurement interval. Fur- 
thermore, the noise introduced by complex measuring equipment can 
obscure the low-level signals encountered in the continuous spectrum* 

The slight discrepancy between the measured and the calculated 
power spectral densities displayed in Figure 8-5 may be largely attributed 
to the non-linearity of the pulse -po sition modulator for the larger modu- 
lating voltages* Consequently^ the standard deviation of the pulse shift 
C is not as large as that determined from the measured standard devia- 
tion of the modulating process* This results in a smaller continuous 
component and corresponding increases in the discrete components as 
compared to the calculated power spectral density, which was based upon 
a linear modulation scheme* 
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Normalized Pulse Width u/T 
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APPENDIX A 



The Pulse-Position Modulator 

A pulse -position modulator was designed and constructed so that 
experimental measurements could be made on the type of PPM waveform 
used for the mathematical model* The design specifications were to 
build a device that sampled a modulating process at a periodic rate and 
generated a pulse train, the pulses of which being delayed in time by an 
amount linearly related to the sampled value* To accomplish this, it 
was decided to build a sampler which had as its "hold" circuit an RC net- 
work with a time constant that resulted in an approximately linear decay 
during the interval under consideration* This linearly decaying voltage 
was then compared to a "zero modulation” reference* When the two 
voltages were equal, a pulse was transmitted to a monostable multivi- 
brator. The mono stable multivibrator was used to give a pulse output of 
the desired shape, and to allow easy variation of pulse width and pulse 
amplitude* See Figure A -2 for pictures taken of the unmodulated and the 
modulated pulse trains at the multivibrator output. 

The circuit diagram is shown in Figure A-l* The variable param- 
eters of the circuit are the pulse width (change Rg), the pulse height 

(change V ), and the sampling rate (change R 1 and pulse rate of the 
cc l 

pulse generator)* The sampling rate is, of course, the average pulse 
repetition rate of the PPM waveform. 



84 



PPM Output 



u 

u 

> 

+ 




tn 



85 



igure A-l. Circuit Diagram of the Pulse-Po sition Modulator. 




Unmodulated Pulse Train 




Pulse train modulated with wideband Gaussian noise. 



Figure A ~2. Performance of the Pulse-Po sition Modulator. 
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